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Abstract. An intriguing potential signature of hydrodynamic behavior in relativistic A+A reactions 
at Relativistic Heavy Ion Collider (RHIC) energies is conical flow induced by fast supersonic 
particles traversing the hot and dense medium. Here I present first results on the evolution of Mach 
shocks in 2 — > 2 covariant transport theory, in a static uniform medium. 
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INTRODUCTION 

3 . 

There has been a lot of recent theoretical interest in conical flow in heavy-ion 
collisionsi! 0, 1, H. 

Important open questions are how the collision dynamics (rapid 
expansion, phase transition, inhomogeneities) and dissipation (viscosity) affect this 
unique flow pattern. 
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Q\ [ It would be natural to investigate these questions using causal dissipative hydrody- 

namics, however, the 3+1D codes required still have to be developed. Here I employ 
a convenient alternative[5], covariant transport theory, which is fully causal and stable, 
can address full 3+ ID, and has a well-defined hydrodynamic limit. 
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MAIN RESULTS 



I consider, as Refs. [@, 0, 0, & 03], Lorentz-covariant on-shell Boltzmann transport 
theory with elastic 2 — » 2 rates, but here couple the system to an external current through 
the transport equation: 

pf^/i =S(x,p 1 )+C[f,f\{x,pi)+C[f at ,f\{x,pi) (1) 

where 

C [frg] = J J (f3g4-flg2)Wn^34 8*(pi+P2-P3-P<l) ■ (2) 
34 

The integrals are shorthands for /■ = / d pi/ (2Ej), the source term S generates the initial 
conditions, while f ext represents the external current. Because the interest is to study the 
theory near its hydrodynamic limit, I take / and f ext to be the phasespace distribution of 
massless "quasi-particles". This ensures that the equation of state e = 3p is close to that 
of the high-temperature plasma in the early stages at RHIC. The transition probability 



W = s(s — 4m 2 )do/dt is adjusted to control the shear viscosity 7] ~ AT / (5o tr ), where 
o tr is the transport cross section. 

In this exploratory study I compute the flow pattern generated by an external "jet" 
moving in the +z direction in a static, uniform, thermal bath (a massless gas of quarks 
and gluons - 3 flavors, 3 colors) with temperature Tq = 0.385 GeV and shear viscosity 
r\ 7a 0. 075.5', where s is the entropy density. The corresponding energy density and mean 
free path are eo ~ 44.7 GeV/fm 3 and A =0.125 fm. Though covariant transport treats 
interactions between jet and medium self-consistently, I here explicitly turn off jet recoil, 
similar in spirit to linear response studies. The jet is created at t = 0. 

In the first 'perturbative' scenario the jet deposits energy and momentum through 
2 — > 2 interactions as encoded in (OQ). This implies dE/dL w dp z /dL at high jet en- 
ergies E ^> T . Typical Debye- screened t— channel processes are quite inefficient at 
energy-momentum transfer to the medium, dE/dL^ (Hjj/hMFpT) \n(ET / . To max- 
imize effects, I therefore take more optimistic isotropic scattering, for which dE/dL ~ 
E/(2Xmfp)- The jet is modeled through a moving sharp sphere profile 

fT-(t,r,z,p) - ^ & (R 2 -( r 2 + (z-t) 2 ))d 2 (p T )d(p z -E) (3) 

of radius R = 0.2 fm, where E is the jet energy, and cylindrical (r,z, 0) coordinates are 
employed. 

Figure \T\ shows the jet-induced change in energy density and momentum density in 
the medium, at t = 2.5 fm, for E = 8.25 GeV and dE jdh « dpJdL « 20 GeV/fm 3 . The 
largest effect is, of course, right at the position of the jet but there clearly is a bow-shaped 
front trailing the jet. In the wake of the jet (z < 2 fm, r < 0.3 fm), momentum flows along 
the jet direction (characteristic "diffusion wake"). These results agree qualitatively with 
ideal hydrodynamic calculations [|4j]. This is quite remarkable because in this calculation 
energy-momentum deposition is not thermal, and there is also a finite, albeit small 
viscosity. 

As an alternative scenario, consider a source that deposits thermalized energy only 
but no momentum {dp z /dL — 0). This can be incorporated through additional thermal 
particle production, i.e., adding to the source term S in (OQ) a contribution 

AS(t,r,Z,p) - v ^l e -[r 2 +(z- Vt ) 2 ]/(2^) e - P /T (4) 

dL 

while putting f ext = 0. Here a Gaussian spatial profile was chosen with a width o = 
0.3 fm, and v is the source velocity. 

Figure [2] shows the change in energy density and momentum density in the medium, 
at t = 2.5 fm, in the "pure energy" scenario for dE/dL ps 75 GeV/fm 3 and v = 0.9c. 
Because the source is turned off for t > 2 fm, the bow-shock is now more pronounced. 
In contrast to the "perturbative" scenario, in the wake region we see momentum flow 
away from the jet. In the regions with the highest momentum density, the flow angle 
is consistent with w 50 degrees expected from the Mach formula cos0 = c s /v (here 
c 2 = 1 /3). These features agree qualitatively with ideal hydrodynamic calculations^. 

For bow shocks to manifest, a low viscosity is crucial. As seen in Fig. [3j the distur- 
bances get largely "washed out" if the shear viscosity is quadrupled to 77 pa 0.3*. 
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FIGURE 1. Distribution of the local energy density (left) and the magnitude of momentum density 
perturbations (right), at t — 2.5 fm for the "perturbative" scenario. Energy density contour lines are at —3, 
— 1, 1, 3, 10 and 40 GeV/fm 3 ; while for the momentum density at 0.5, 1, 2, 3, 4, and 20 GeV/fm 3 . Arrows 
indicate the local momentum density vectors. 




For the modest number of test particles and events in this exploratory study, there 
were no clear signals in the azimuthal distribution in the x — z plane, dN /dydfy, and the 
angular distribution relative to the jet axis dN /d cos 6 (even when weighted by powers 
of momentum). Ideal hydrodynamic calculations also find very small effects, confined 
to thermal tails p/T ~ 20 Q4£| . It remains to be seen whether the dynamics and self- 
consistent jet-medium coupling can generate appreciable signals in heavy-ion collisions. 
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FIGURE 3. Distribution of the local energy density perturbations for the "perturbative" (left) and "pure 
energy" scenarios (right), same as Figs.Qlleft) and Figs.|2jleft), but with four times larger mean free path. 



CONCLUSIONS 

In this work I investigate Mach shocks in 2 — > 2 covariant transport theory, in a static 
uniform medium. If the shear viscosity to entropy density ratio is very low, r\/s~ 0.075, 
the results are in qualitative agreement with ideal hydrodynamic calculations [4]. This 
demonstrates the feasibility of utilizing covariant transport in future 3+ ID conical flow 
studies to incorporate viscosity and a self-consistent coupling between jet and medium. 
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